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ABSTRACT 

We present spectral line observations of the ground state transitions of hydroxyl(OH) toward 
supernova remnant (SNR) IC443 carried out with the Green Bank Telescope. Weak, extended 
OH(1720 MHz) maser emission with OH(1667, 1665, 1612 MHz) absorption is detected along the 
southern extent of the remnant where no bright compact maser sources have been previously 
observed. These newly detected SNR-type masers are coincident with well-known molecular 
clumps and a ridge of shocked H2 emission indicative of the SNR shock front interacting with 
the adjacent molecular cloud. Simultaneous observation of all four ground-state transitions of 
OH permits us to fit physical conditions of the shocked gas at the interaction site. A simple two- 
component model for the line profiles yields the physical parameters for detected regions of maser 
emission including excitation temperature, OH column density and filling factor. Observed line 
profiles suggest the shock is largely propagating toward the line-of-sight in the region of these 
newly identified weak masers. The implications of shock geometry and physical parameters 
in producing spatially extended OH maser emission in SNRs are explored. We also present 
VLA radio continuum observations at 330 MHz for comparison with OH line observations of the 
remnant. 

Subject headings: ISM: general — shock waves — supernova remnants: individual(IC443) — 
Facilities: GBT, VLA 



1. Introduction 



The supernova remnant(SNR) IC443 (G189. 2-1-3.0) is one of the best sources in which to study the 
physical processes of a supernova remnant interacting with adjacent molecular clouds. Located at a distance 
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of 1.5 kpc in a direction largely free of the confusion present toward the Galactic plane, it was well suited 
for some of the earliest studies of molecular tracers of shock excitation (DeNoyer 1979; Huang, Dickman & 
Snell 1986; Burton 1987; Tauber et al. 1994). Evidence of IC443's interaction was first observed as large 
line- width, negative velocity neutral hydrogen(HI) by DeNoyer (1978). Subsequent carbon monoxide(CO) 
observations showed a slow shock had been induced across a ridge of molcciilar material in the southern half 
of the remnant (DeNoyer 1979). Notably, the ratio of OH/CO column densities was observed to be enhanced 
by two orders of magnitude compared with preshock cloud. Huang, Dickman & Snell (1986) identified a 
series of shocked CO clumps (designated A-H) along this southern ridge, with sizes ranging from 0.5 to 5 
pc. Further high-resolution study revealed shocked molecular line emission tracing out a ring of perturbed 
molecular gas about the center of the remnant (Dickman ct al. 1992; van Dishocck, Janscn & Phillips 1993). 

As a rich, well-studied laboratory for shock chemistry, IC443 was again a prominent target for the study 
of OH following the realization by Frail, Goss & Slysh (1994) that maser emission from the 1720.530 MHz 
transition is a signpost of interaction between SNRs and molecular clouds. Theoretical considerations yield 
strict constraints on the conditions necessary for SNR OH(1720 MHz) masers: moderate temperatures of 
50-125 K, densities of 10^ cm""^, and OH column densities of the order 10^^-10^'' cm~^ (Lockett, Gauthier 
& Elitzur 1999). Masers are thought to be coUisionally excited when a non-dissociative shock propagates 
into a sufficiently dense molecular cloud. 

While these conditions are satisfied throughout much of the shocked ring of molecular material present 

in IC443, only a single molecular clump G near the projected center of the remnant hosts identified masers 
(Frail et al. 1996). Subsequent high-resolution study has used this species of maser as a unique probe 
of post-shock magnetic field, as well as density and temperature at arcsecond resolution (Claussen et al. 
1997). Hoffman ct al. (2003) has resolved the brightest compact masers in clump G, revealing 60 AU cores 
surrounded by diffuse halos extending up to 300 AU in size. While the entire shocked clump appears capable 
of producing maser emission only a small fraction does (~10%). However, such studies of IC443 remain 
limited to clump G. 

The lack of observed compact OH(1720 MHz) emission toward the southern ridge of IC443 underscores 
the importance of geometry in the amplification of OH(1720 MHz) emission. Claussen ct al. (1997) suggests 
that only when the shock is transverse to the line of sight can OH(1720 MHz) masers be observed; it is 
in this geometry that maximal velocity coherence and pathlength occur. As the shock propagates into the 
adjacent cloud, water molecules produced by the shock are dissociated into OH indirectly by thermal X-rays 
emitted from hot gas in the SNR interior (Wardle 1999). This OH enhancement (OH/H2>10-'^) behind 
the shock can be detected as thermal absorption against the bright background continuum of the remnant, 
regardless of shock geometry. Absorption profiles show the velocity dispersion of the gas toward the line of 
sight, which permits shock geometry traced by OH to be determined. 

Motivated by this unique opportunity to probe very specific conditions just behind the shock front we 
have observed the four ground state transitions of OH over the entire extent of IC443. Sections 3.1 and 
3.2 present our detections of OH(1720 MHz) maser emission and ground-state OH absorption, respectively, 
across the southern portion of the remnant. Simple modeling of the identified maser sources is given in 
§3.3. Kinematics of the masing regions and and the nature of extended maser emission are discussed in §4 
in relation to the global shock interaction picture. Our 330 MHz radio continuum image of IC443 is used to 
place the shock interaction picture presented by line observations in context with the radio morphology of 
the remnant. 
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2. 



Observations 



2.1. Spectral Line Data 



The Green Bank Telescope (GBT) of the National Radio Astronomy Observatory^ was used on 2003 
August 15 and October 7 to observe the four ground-state transitions of OH at 1612.231, 1665.4018, 1667.359 
and 1720.530 MHz in left- and right-circular polarizations. The observations presented here are part of a 

larger observing program which mapped all four ground-state transitions of OH for 21 remnants. A more 
detailed account of the results of this survey will be given in a future publication. 

A map of the entirety of IC443 was constructed using a grid of 210 pointings which extends 51'x48' 
about the projected center a,S{J2000) = 6''16™56% -F22°29'8". The GBT has a beam-size of 7.'2 at 1.7 GHz, 
so pointings are spaced by 3.'3 to correspond to the maximal spatial resolution for Nyquist sampling. Each 

pointed spectra has a RMS noise level of 15-20 m.Jy channel" -'^ with a velocity resolution of 0.28 km s~^ over 
the 12.5 MHz bandwidth. Calibration via in-band frequency switching limits the effective velocity coverage 
to 400 km s~^ (1.17 MHz) about the rest frequency of each line. 

For modeling it was necessary to obtain the continuum level from our observations. The difference of 
the system temperature measured on source compared to our blank sky position {AT gys) is used to derive 
the continuum temperature (Ts^c) using the equation: 



where rjB is beam efficiency (0.94), r is atmospheric opacity (0.011 at L-band), and AA = l/sin(Osrc)~ 
l/sin(Osfej,) where O is the elevation of observations on source and on blank sky. We estimate the measured 
continuum temperature to be accurate within 10%. Calibration, image processing and profile fitting utilized 
the AIPS-I— |-^ and MIRIAD^ astronomical processing systems. 



The 330 MHz radio continuum observations of IC443 were obtained as part of an unrelated project to 
survey 7-ray sources from the EGRET catalog (Hartman et al. 1999 and references therein). One of these 
sources, 2EG J0168-I-2234, lies close enough to IC443 with respect to the VLA's primary beam (~2?5 at 
330 MHz) so as to provide a good serendipitous image. The observations were obtained in March and July 
of 1996 using the VLA's C and D configurations, respectively. The data were acquired in 2 IFs centered at 
327.5 and 333 MHz in each of two circular polarizations. Since the observations were part of a survey, the 
on-source integration time was limited to approximately 40 and 30 minutes in the C and D configurations, 
respectively. While limited, these cycling-snapshot multi-configuration observations nonetheless provided a 
sufficient range of spatial frequencies to recover most of the structure in this large (~45') object. 

Post-processing of low-frequency VLA data uses procedures similar to those at higher frequencies al- 
though certain details differ. Phase calibration serves to compensate for ionospheric phase effects but because 
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■^MIRIAD is the Mutlichannel Image Reconstruction, Image Analysis and Display system of the Berkeley-IUinois-Maryland 
Association (BIMA) 



ATsys = 3AA + r]BTsrcexp{-TAsrc) 



(1) 



2.2. 



Continuum Data 
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the ionosphere is fairly coherent across the VLA in these configurations single observations of 0137+331 
(3C48) and 0542+498 (3C147) were sufficient for both flux density, bandpass and initial phase calibration. 
In order to combat RFI, the data were acquired in multi-channel continuum mode (32 channels per polariza- 
tion) at higher spectral resolution than needed for imaging. Furthermore, the large field of view requires full 
three-dimensional dcconvolution to remove the effects of the VLA's non-coplanar baselines. To achieve this, 
the data were edited and calibrated in standard AIPS^, while NRAO's Software Development Environment 
polyhedron algorithm "DRAGON" (Cornwell & Perley 1992) was utilized for performing successive loops 
of self-calibration and imaging. Comparison of the position of NVSS sources detected in our image with 
their positions in the NVSS catalog (the latter good to ^1") showed that the astrometry of the radio map 
is accurate to at least a few arc-seconds. 

The final 330 MHz radio continuum image, uncorrected for primary beam attenuation, is shown in 
Figure 5. The synthesized beam is 74"x64" at a position angle of 48 degrees, and the rms noise is 6 mJy 
beam^^. After removal of unrelated point sources, the peak brightness of the SNR is 870 mJy beam~^ and 

the total integrated flux is 182 .Jy after primary beam correction. Because cycling snapshots are notoriously 
limited in their ability to recover emission from such an extended object these flux measurements should be 
taken as a lower limit. Comparison with single dish measurements of this object compiled by Erickson & 
Mahoncy (1985) give flux densities well over 200 Jy, even at somewhat higher frequencies. However, our 330 
MHz continuum image does reproduce most of the structure revealed in the deeper 327 MHz WSRT map 
(Braun & Strom 1986) and is more than adequate to provide a useful finder atlas for referencing our higher 
frequency OH observations. 



3. Results 

3.1. OH(1720 MHz) Maser Emission 

We have identified three sites of OH(1720 MHz) emission in IC443 with ^1 km s~^ line- widths indicative 
of stimulated emission. At the T.2 resolution of these observations we identify a bright (Sp=3.6 Jy) maser 
feature associated with molecular clump G, and two spatially extended maser features of Sp«0.15 Jy which 
are coincident with molecular clumps B and D. Table 1 gives the Gaussian fits to these emission features. 
Here our maser designations follow those of molecular clumps presented in Huang, Dickman & Snell (1986, 
see Figure 1). 

The bright maser G has been previously detected with 1.5" resolution VLA observations as six compact 
maser spots all with velocities near -4.5 km s~^ (Claussen et al. 1997). All six compact masers fall well 
within one GBT pointing. The single dish measured flux density of 3.64 Jy is consistent with the VLA 
measured flux density of 3.89 Jy, accounting for single-dish flux calibration errors of 10% . The implication 
is that all maser emission from this clump is present on the very small scales resolved by interferometric 
observations. Right circular polarization is detected at a level of 13%, consistent with the earlier results of 
Claussen et al. (1997) and Hoffman et al. (2003). 

For the first time OH(1720 MHz) maser emission is detected along the southern extent of IC443 at a 
level 20 times weaker than the bright compact maser G. Clumps B and D are observed with a peak brightness 
of 158 and 145 mJy and at a velocity of -6.14 and -6.85 km s~^, respectively. Claussen et al. (1997) were 
not able to detect weak maser features toward these clumps down to 5 mJy beam~^ channel"^ with high 
resolution VLA observations. Maps of the velocity integrated line fiux density (moment zero) in Figure 2 
show the newly detected emission features extending several arcminutes along the southern ridge of IC443 
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following the nonthermal radio continuum at 330 MHz. Circular polarization was not detected for these 

weak, extended features. 

We consider two possible explanations for the weak 1720 MHz emission along the southern ridge. First, 
the observed emission could be due to multiple compact masers coincident with clumps B and D. However, 
previous VLA observations would have been sufficient to detect these sources if the masing regions are smaller 
than 1' in size, as was observed toward clump G. Therefore we suggest that these newly identified masing 
regions are spatially extended, having properties similar to the weak, extended regions of maser emission 
observed in remnants W28, 3C 391, G359. 1-0.5, G357.7+0.3 and G357.7-0.1 (Yusef-Zadeh, Uchida & Roberts 
1995; Yusef-Zadeh et al. 1999). The narrow line-widths observed (<2 km s~^) are indicative of non-thermal 
maser emission. Brightness temperatures for these extended masers are ^7, K, uncharacteristically low 
for compact masers (Tb «10'' K) but consistent with the several arcminute scales and lower brightness 
temperatures (Ts »10-2500 K) of extended maser emission identified by Yusef-Zadeh et al. (1999). The 
implications of extended maser emission are discussed in §4.2. 

It is unlikely that the newly identified maser emission is a result of time-variable masing that has 
brightened significantly since VLA observations were conducted. Goss (1968) first detected OH(1720 MHz) 
maser emission in W28 and no time- variability has yet been observed for that source (Frail et al. 1994). We 
observe no change in the brightness of emission from clump G compared to previous observations over an 8 
year period (Hoffman et al. 2003; Claussen et al. 1997). 

3.2. Thermal OH(1667,1665,1612 MHz) Absorption 

Simultaneous observations of all four transitions of OH confirm that masing occurs only at 1720 MHz and 
traces the kinematic structure against the nonthermal continuum of the remnant. The strongest absorption 
and clearest velocity features are observed in the 1667 MHz transition. Figure 3 shows line profiles of the 
four observed transitions toward maser clumps G, B and D. The absorption is found strongest in the 1667 
MHz transition, and weakest at 1612 MHz. No significant absorption is detected in the 1720 MHz transition. 

The absorption profile toward clump G can be well-fit by two Gaussian components; (1) a stronger, 
narrow line with AV«5 km s^^ at ^ hSR~~'i km s^^; and (2) a weaker, wide line with AV«20 km s~^ at 
km s~^. This profile is characteristic of transverse acceleration, where the absorption is strongly 
peaked near the systemic velocity of the remnant, and is somewhat broadened by thermal excitation of 
OH gas from the passage of the shock (Figure 3a). Furthermore, there is no clear delineation in velocity 
between pre- and post-shock absorption, though close inspection of the spectra of clump G reveals a slight 
asymmetrical broadening toward negative velocities indicating part of the shock is inclined slightly toward the 
line of sight. Other molecular species observed toward clump G show the same characteristic profile (Burton 
1987; Dickman et al. 1992; van Dishoeck, Jansen & Phillips 1993; Tauber et al. 1994). A transversely 
propagating shock is also needed to maximize the coherent path-length and amplification of the observed 
maser (Claussen et al. 1997). 

Absorption profiles throughout the southern extent of the remnant show broad, asymmetrically blue- 
shifted line- widths of 30-40 km s~^. Figures 3b and 3c show line profiles toward masers B and D, respectively. 
Strongly asymmetric blue-shifted broadening of the post-shock absorption profile suggests that the shock is 
accelerating the molecular gas toward the observer. In this geometry gas in all phases of shock acceleration is 
present along the line of sight causing this spectral delineation of velocity components relative to the systemic 
velocity of the remnant. A narrow line at Vlsr ~ -3 km s~^with a width AV < 5 km s~^delineates the 
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ambient unshocked gas. The shocked gas along the southern ridge is fit by two Guassian components: (1) 
a strong, broad Hne with AV ~ 10 km s~^at Vlsr ~ -10 to -12 km s~^; and (2) a shallow, very broad 
line with AV « 20 to 30 km at Vlsr ~ -15 to -25 km s~^. This decomposition of the asymmetric 
broadening into two Gaussian components can be non-unique, but is offered here as a means of first-order 
insight into shock dynamics. 

Thus, the shock geometry is inferred to be transverse to the line of sight through clump G, and toward 
the line of sight along the southern ridge that includes clumps B and D, where extended masers are detected. 
The association between shock geometry and extended maser kinematics is discussed further in §4.1. 

There is good agreement between the profiles of OH and other molecular shock tracers, particularly CO, 
observed at comparable resolution. Snell, et al. (2005) observed H2O, CI, 1^00(5-4), 12CO(1-0), and HCO+ 
for clumps B and G (see Figures 1 and 3, respectively in that work). -^^00(5-4), thought to be created by 
the shock, mirrors the kinematic picture shown by OH absorption in clumps B and G. However, for both 
clumps H2O profiles extend beyond the largest velocities at which we detect OH absorption: in clump G 
H2O emission is is observed out to -30 km s~^and in clump D out to velocities of -80 km s~^. Snell et al. 
(2005) have suggested both a fast J-typc and slow C-typc shock are needed to explain the observed shocked 
molecular species along the southern extent of IC443. The observed broader velocity components only for 
H2O is consistent with some H2O being liberated behind a J- type shock, while OH is mainly created behind 
the C-type shock front. It is interesting to note that given the observed OH column densities (Table 2) 
the column of H2O is orders of magnitude less than predicted by theory. This apparent lack of H2O is not 
accounted for in any present theory for the production of OH(1720 MHz) masers. 



In this section we estimate the parameters of the gas responsible for producing the observed line profiles 
using a simple two-component model. As shown in Figure 4, our model fitting includes a maser component 
and a component to fit the broad absorption arising from the shocked gas. We do not include ambient 
gas absorption in these simple models. Each component is characterized by its line-of-sight OH column 
density A^oh, kinetic temperature T^, mean and FWHM velocities, vq and Av, and a beam filling factor /. 
The H2 density is fixed at 10^ cm^"^, typical of what is expected in OH(1720 MHz) mascr-cmitting regions 
(Lockett, Gauthier & Elitzur 1999). Excitation by collisions and partial trapping of photons then determines 
the excitation temperature Tex of the 4 transitions and the line-center optical depth. Model fits are done 
independently of the Gaussian fits used for the discussion of kinematics in earlier sections. 

The observed intensity for a particular line is 



where Iqv is the background continuum intensity and B^, is the Planck function. Here the obscuration of 
one component by another is neglected. Assuming that |Tex| ^ 0.08 K, the Rayleigh-Jeans approximation 
holds, and eq. (2) can be recast as 



3.3. Emission and Absorption models 



= (1 - /)/o. + / [lo. exp(-T,) + B,(Tex)(l - cxp(-T,))] 



(2) 



T{v) - To = /(Tex - To)(l - exp(-T(^;))) 



(3) 



We write the optical depth as 




(4) 
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where tq is the optical depth at line center and a = At;/(81n2)^/^ and we have adopted a Gaussian line 
profile. 



The optical depth at line center is 



To 



ln2V/^ X^ANoiig2Xi 



4n Av 



91 



1 — exp 



AE 



(5) 



where gi,xi and g2,X2 arc the degeneracies and fractional populations of the lower and upper states, and 
A, A and AE are the wavelength, Einstein A value, and energy difference for the transition. 

For the four levels in the rotational ground state of OH, usually |Tcx| S> AE/k ^ 0.08 K; in addition 
the higher rotational states are not significantly populated. This means that the term in the square brackets 
in eq. (5) can be approximated as AE/kTex and also that the four levels in the ground rotational state are 
very nearly populated according to their degeneracies (ie. xi^2 ~ 51,2/16). Thus 



To 



A^,u(10''c:ui-2) 
' Tex(K)At;(km/s) 



(6) 



where the constant a is 0.4540, 2.3452, 4.2261, and 0.4846 for the 1612, 1665, 1667 and 1720 MHz lines 
respectively. Note that Tex and tq are negative for an inverted transition. 

Eqs (3), (4) and (6) form the basis of a simple model for the observed lines, with parameters Aqhi ^fe, 
/, Vq, Av. The excitation temperatures are determined by an LVG model of the coUisional excitation of 
the OH molecules within a clump, including the effect of line overlap and radiative transfer in the maser 
line, following Lockett, Gauthier & Elitzur (1999). Given the uncertainty inherent in the single-dish data, 
we do not perform a detailed fit to the line profiles (which in any case arc not Gaussian); instead we have 
adjusted parameters until a reasonable match is found. The derived parameters should therefore be regarded 
as illustrative rather than definitive. 

To match the observed profiles two Gaussian components are required: a narrow component to match 
bright 1720 MHz maser emission, plus a broad component responsible for the bulk of absorption in the 

other lines. Fits of our simple models for clumps B, D and G are presented in Figures 4-a,b,c respectively. 
The net contribution of components is shown as a bold dark line against the observed profile. The broad 
component has a peak optical depth of order unity for Clump G and appears optically thin for Clumps B 
and D (see Table 2). For the line of sight to clumps B and D a single broad component with parameters very 
similar to the broad component for the clump G model, except with a significantly lower column density, 
successfully reproduces the line profiles in all four transitions. Given the observed OH(1720 MHz) narrow 
emission component, the column densities derived from our simple models are as expected from models of 
OH production in shock waves (Wardle 1999) and collisional excitation of OH masers (Lockett, Gauthier & 
Elitzur 1999). 

For the narrow component the main lines do not appear in the expected rough 5:9 ratio for optically 
thin absorption with similar excitation temperatures. We observe Tex ^ Tq w iOK, and the brightness 
temperature at line center only drops by 0.3 K. The relatively bright and narrow 1720 MHz maser component 
requires a much smaller filling factor to avoid introducing sharp absorption feature into the absorption lines. 
We therefore conclude the line is optically thick with a beam filling factor / ~ 0.01 for this component. We 
note that some confusion between pre- and post-shock gas may be unavoidable for modeling of clump G 
given the overlapping velocities of ambient cloud and post-shock gas. For clumps B and D the separation 
between pre- and post-shock velocity components is sufficient to avoid significant confusion between the two. 
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This is apparent from the significant absorption present at -3 km s which extends below the dark Hne 
indicating the simple model fit in Figures 3b and c. 

3.4. OH Abundances 

To further demonstrate that OH is enhanced behind the passage of the SNR shock, as predicted by 
theory we can compare pre- and post-shock columns and abundances. The column density of OH can be 
estimated using 

TVoH = 2.2785 X 10" T^^ j Ty dv cm'^ (7) 

for the 1667 MHz line (Crutcher 1977). For clump B, in which the separation between pre- and post-shock 
components is clearest, we find Nqh = 1.2x10^^ Tex cm~^ for the pre-shock gas. Here we have used the 
Gaussian fit to the pre-shock component presented in §3.2 to derive Tpeaft=0.026 and AV=2 km s~^. From 
our simple modeling of the broad post-shock gas we find Nqh = 1.1 x 10^^ cm~^ and Tex = 7.5 K at 1667 MHz. 
Assuming the excitation temperature for the pre-shock gas is no more than that found for the shocked gas, 
we find an enhancement in OH column greater than two orders of magnitude. This enhancement in column 
density is indicative of the production of OH behind the shock front given that the shock through clump B 
is propogating mostly along the line-of-sight and both the pre-shock and broad component post-shock gas 
are thought to fill the beam. 

To calculate OH abundances we used characteristics of the pre- and post-shock gas derived for clump 
B by van Dishoeck, Jansen & Phillips (1993). The pre-shock column of H2 is found to be 2x10^^ cm~^ 
yeilding OH/H2=4.5x 10~*. For the post-shock gas the CO column density is found to be l.lxlO^'' cm~^ for 
a gas with T=80 K and niif2=10^ cm~^ and AV = 20 km s~^, comparable to the broad component we have 
modeled for clump B. Assuming CO/H2 is 10"'', typical of interstellar clouds, the post-shock abundance of 
OH is 6x 10~^ indicative of an OH enhancement exceeding two orders of magnitude. We note that molecular 
clumps in IC443 are typically ~1' in extent (van Dishoeck, Jansen & Phillips 1993), which is significantly 
smaller than the 7'.2 beam of these observations. Beam dilution may be significant, up to a factor of 50, 
however we argue that such a large dilution value is unlikely. Large-scale CO observations show the pre-shock 
gas is extended over regions greater than the GBT beam (Seta et al. 1998). Our simple modeling of the 
post-shock gas finds a filling factor of order unity for the broad component of OH absorption, consistent 
with multiple shocked clumps of arcmimite size present within a single pointing. Even accounting for the 
uncertain effects of beam dilution, it is clear that OH is enhanced by the passage of the SNR shock. 

4. Discussion 
4.1. Kinematics 

OH(1720 MHz) maser amplification favors a very specific geometry for the SNR-MC interaction in which 
the shock propagates across the sky, transverse to the line of sight. In this geometry the velocity gradient 
of the masing gas is minimized and the largest coherent path lengths are achieved (Frail et al. 1996). Thus 
the observed line-of-sight velocities of SNR- type masers are commonly found within a few km of the 
systemic velocity even when distributed throughout the remnant. 

The masers associated with 1C443 all lie within 2 km s~^ of the systemic velocity (-5 km s~^), and 
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the compact maser emission of clump G has been associated with a shock propagating transverse to the 
hne of sight as is expected (Claussen et al. 1997). OH absorption profiles toward clumps B and D show a 
clear delineation between pre- and post-shock velocity components with asymmetrically blue-shifted shock 
broadening, indicative of shock acceleration toward the line of sight. This face-on geometry is observed along 
the entirety of the southern ridge of IC443 tracing the shock interface with adjacent molecular material. In 
this face-on geometry the OH pathlength and maser intensity are greatly reduced, and masing gas is present 
at distinctly different velocity than the bulk of the gas observed in absorption. 

To reconcile these disparate kinematic results, it is necessary to consider the global kinematic structure 
revealed by high resolution molecular line observations. Dickman et al. (1992) mapped the entirety of IC443 
tracing shocked gas in the J=l-0 lines of CO and HCO"*". The observed shocked clumps arc found in a ring 
~9 pc across, tilted by ~50° from the line of sight. Observations of ^^CO 3-2 transition by van Dishocck, 
Jansen & Phillips (1993) resolve molecular clumps A-H into smaller sub-clumps embedded in a sinuous ring 
of molecular material. At 1.6 GHz the GET beam samples regions with 7' resolution, and may include 
multiple molecular sub-clumps within one pointing. 

For all positions in which 1720 MHz maser emission is observed, comparison with CO 3-2 spectra 
indicates the existence of transverse shocks in which an optimal geometry for OH(1720 MHz) maser emission 
can arise. For a comparison of CO 3-2 with our OH observations we note that Figure 3 of van Dishoeck, 
Jansen & Phillips (1993) presents spectra of CO sub-clumps RlO-13 contained within clump B, R17-19 
within clump D, and Rl-2 within clump G. For mascrs B and D, while the brightest sub-clumps show highly 
blue-shifted profiles, weaker sub-clumps are also identified within each GBT pointing that show CO 3-2 
profiles indicative of transversely directed shocks. At yet higher resolution, CO 3-2 sub-clumps within clump 
B show drastic variations in observed line profiles on scales less than 1' (Figure 5 of van Dishoeck, Jansen 
& Phillips 1993). Small scale variations in shock stucture are consistent with profile modeling in §3.3. Both 
a broad and narrow component are needed to explain the observed profiles, representing the face-on and 
transversely accelerated clumps, respectively. Modeling of the masing gas also indicates OH column densities 
are a factor of 5 to 20 lower for masers B and D than for maser G. Maser emission arises in the low-column 
density, transversely shocked clumps throughout the southern ridge. 

Globally, blue-shifted OH absorption profiles are never foimd to exceed « -20 to -30 km s^^. This is 
consistent with the ring of molecular material undergoing a non-dissociative C-type shock. As shown by 
Lockett, Gauthier & Elitzur (1999) and Wardle (1999), only C-type shocks are capable of producing the 
conditions necessary for the production of OH(1720 MHz) masers. This is in contrast with observed 0-H2O 
profiles which extend out to -80 km s~^for clump B. It has been suggested that a faster J-type shock is 
needed to explain these large velocity offsets for shocked molecularl material and the observed [C II] emission 
(van Dishoeck, Jansen & Phillips 1993; Snell et al. 2005). However, OH appears to only be present at lower 
velocities consistent with a slow shock. 



4.2. Extended Maser Emission 

As discussed in §3.1, the newly detected masers B and D are thought to be spatially extended regions 
of masing gas. In previous cases compact maser sources have been identified with surrounding extended 
halos of maser-like emission at significantly lower brightness temperatures: typically Tb~1-2500 K whereas 
compact maser sources are typically Tb«5x10* K (Yusef-Zadeh et al. 1999). Such a low brightness level 
and a large spatial extent would have made previous interferometric observations with the VLA insensitive 
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to such extended maser emission. 

Gaussian fits to these newly detected features show spatial extension over a few GBT beam widths. 
Those large scale masers suggest a low-gain population inversion of OH is created over extended regions of 
the ambient cloud by the passage of the shock front. This is supported by the correlation between regions 
of extended maser emission and shocked H2 emission indicative of the SNR shock front interacting with 
the adjacent molecular cloud. Such is the case for the newly identified extended maser emission in IC443, 
which is present along the southern ridge of H2 emission (Burton 1987) As discussed in §4.1, the molecular 
environment at these interaction sites is clumpy on scales smaller than the resolution of our observations 
(Tauber et al. 1994), so further observations which can resolve the weak, extended regions of masing are 
needed to probe the shock-cloud interaction on scales comparable to the dumpiness of the medium. 

The newly detected extended masers in IC443 are the first such instance of weak, extended maser 
emission with no associated compact masers. A velocity gradient of a few km s^^ is present over the expanse 
of the extended maser emission suggesting that for large-scale masers acceleration is unlikely to be entirely 
perpendicular to the line of sight, as is thought to be the case for compact masers (Lockett, Gauthier & 
Elitzur 1999). IC443 may not be unique in this respect; extended masers in 3C 391 appear at significantly 
different velocities from compact masers. Further modeling of maser amplification under different shock 
geometries is necessary to better understand the stimulation of these extended masers. 

5. Conclusions 

Sensitive observations of the four ground-state transitions of OH over the entirety of IC443 lead to the 

following conclusions. 

1. Thermal OH absorption is detectable against the background continuum throughout the southern 
extent of 1C443, and reveals the direction of shock propagation. The western extent (1C443G) is undergoing 
transverse motions while bulk motions throughout the southern extent (IC443B,D) indicate a line of sight 
acceleration. 

2. Extended maser emission is detected along the southern ridge of IC443 toward molecular clumps 
B and D with low brightness temperatures and velocities systemic to the remnant and associated ambient 
medium. The disparity between observed maser emission and absorption profiles is indicative of cloud 
dumpiness and kinematic variations on scales less than 7.'2, the resolution of these observations. 

3. The observed OH line profiles cannot be fit with a single shock component. Two component modeling 
derives vq, Av, Noh, T;;, / and T^x for detected regions of maser emission, consistent with theoretical 
constraints for OH(1720 MHz) masers in supernova remnants. 

Further investigation of weak extended masing must be conducted at higher resolutions to constrain the 
spatial scales and shock geometries in which extended maser emission occurs. 

We are very grateful to Ron Maddalena and Jim Braatz for their assistance with observations at Green 
Bank and helpful discussions. Support for this work was provided by the NSF through awards AST-0307423, 
and GSSP 2-0005 from the NRAO. Basic research in radio astronomy at the Naval Research Laboratory is 
supported by the Office of Naval Research. 



- 11 - 



REFERENCES 

Braun, R. & Strom, R.G. 1996, A&A, 164, 193 
Burton, M. 1987, Q. Jl R. Astr. Soc. 28, 269 

Claussen, M.J., Frail, D.A., Goss, W.M., & Gaume, R.A. 1997, ApJ, 489, 143 

Cornwell, T.J. & Perley, R.A. 1992, A&A, 261, 353 

Crutcher, R.M. 1977, ApJ, 216, 308 

DeNoyer, L.K. 1978, MNRAS, 183, 187 

DeNoyer, L.K. 1979, ApJ, 232, L165 

DeNoyer, L.K. 1981, ApJ, 246, 37 

Dickman, R.L., Snell, R.L., Ziurys, L.M. & Huang, Y. 1992, ApJ, 400, 203 
Erickson, W.C. & Mahoney, M.J. 1985, ApJ, 290, 596. 
Frail, D.A., Goss, W.M., & Slysh, V.L 1994, ApJ, 424, LI 11 

Frail, D.A., Goss, W.M., Reynoso, E.M., Giacani, E.B., Green, A.J. & Otrupcek, R. 1996, AJ, 111, 1651 
Goss, W.M. 1968, ApJS, 15, 131 
Hartman, R.C. et al.. 1999, ApJS, 123, 79 

Hoffman, I.M., Goss, W.M., Brogan, C.L., Claussen, M.J. & Richards, A.M.S. 2003, ApJ, 583, 272 
Huang, Y.L., Dickman, R.L. & Snell, R.L. 1986, ApJ, 302, L63 
Lockett, R, Gauthier, E. and Elitzur, M. 1999 ApJ, 511, 235 

Seta, M., Hasegawa, T., Dame, T.M., Sakamoto, S., Oka, T., Handa, T., Hayashi, M., Morino, J., Sorai, K. 
& Usuda, K.S. 1998, ApJ, 505, 286 

Snell, R.L., HoUcnbacli, D., Howe, J.E., Neufeld, D.A., Kaufman, M.J., Melnick, G.J., Bergin, E.A. and 
Wang, Z. 2005, ApJ, 620, 758 

Tauber, J.A., Snell, R.L., Dickman, R.L. and Ziurys, L.M. 1994, ApJ, 421, 570 

van Dishoeck, E.F, Jansen, D.J. & Phillips, T.G. 1993, A&A, 279, 541 

Wardle, M. 1999, ApJ, 525, LlOl 

Yusef-Zadeh, F., Uchida, K.I. & Roberts, D.A. 1995, Science, 270, 1801 

Yusef-Zadeh, F., Goss, W.M., Roberts, D.A., Robinson, B. and Frail, D.A. 1999, ApJ, 527, 172 



This preprint was prepared with the AAS lATJiX macros v5.2. 



- 12 - 



Table 1: OH(1720 MHz) Emission Gaussian Fits 
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o / // 
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O 


B 


06 16 40.3 


+22 23 06 


158(22) 
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G 


06 16 47.0 


+22 32 01 
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Errors to fits are shown in parenthesis. 



Table 2: Model parameters 
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0.27 
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-0.39 
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Fig. 1. — Gray-scale continuum image of IC443 at 330 MHz between -300 and 700 mJy beam ^ at a 
resolution of 74"x64". Positions of maser emission are demarcated with crosses which extend one GBT 
beam-width {7'.2). 
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Fig. 2. — Top left: 330MHz continuum map of IC443; Top right: the velocity integrated line flux density 
(moment zero) map of maser G at -4.55 km s^^. Bottom right: moment zero map of maser B at -6.14 km 
s^^; Bottom left: moment zero map of maser D at -6.85 km s~^. Contours of 330 MHz continuum are 
present for all figures at 0.01, 0.05, 0.1 and 0.3 times the peak continuum level of 177 mJy beam^^. Note 
that due to the brightness of maser G its emission wings are seen in all three moment zero maps. 
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Fig. 3. Spectra of all four ground-state transitions of OH from pointings taken toward three clumps with 
detected 1720 MHz emission, (a) Clump G at a,(5(J2000) = 6^16'"45!4, +22°32'16". A strong, narrow 
emission feature is seen at 1720 MHz, whereas broad absorption is seen at the other transitions of OH, 
indicating a transverse shock, (b) Clump B at a,^(J2000) = 6''17'"29!4, +22°23'38". A narrow emission 
feature is seen at 1720 MHz just above the 3<t level, whereas broad absorption is seen at the other transitions 
of OH, with a sharp absorption line representing the preshock gas at Vlsr = -3.21 km s~^. (c) Clump 
D at a,^(J2000) = 6''17"'57!6, +22°24'34". A narrow emission feature is seen at 1720 MHz just above 
the 3(7 level, whereas broad absorption is seen at the other transitions of OH, with a sharp absorption line 
representing the pre-schock gas at Ylsr = -4.32 km s~^. 
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Fig. 4. — Spectra of all four ground-state transitions of OH from pointings taken toward three clumps 
with detected 1720 MHz emission. The solid line shows the best fit model for the observed spectrum, (a) 
Clump G at a,5(J2000) = 6^16™45?4, -|-22°32'16". A strong, narrow emission feature is seen at 1720 MHz, 
whereas broad absorption is seen at the other transitions of OH, indicating a transverse shock, (b) Clump 
B at a,(5(J2000) = 6'^17'^29!4, +22°23'38". A narrow emission feature is seen at 1720 MHz just above 
the 3(7 level, whereas broad absorption is seen at the other transitions of OH, with a sharp absorption 
line representing the preshock gas at Vlsr = -3.21 km s~^. (c) Clump D at a, i5(J2000) = 6'^17™57;6, 
+22°24'34". A narrow emission feature is seen at 1720 MHz just above the 3(t level, whereas broad absorption 
is seen at the other transitions of OH, with a sharp absorption line representing the pre-schock gas at Ylsr 
= -4.32 kms-^ 



